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1 Introduction 



We want to report in this article the calculation of transition amplitudes of multi-magnon 
scattering in the ferromagnetic Heisenberg XXX-chain with an inhomogeneous distribution of 
the coupling constant. 

The Hamiltonian of the model under consideration is given by 

H = H hom + H inh (1) 
J N 

H hom = 7 X [Wn+l ~ 1] (la) 
4 n =-N 

1 N 

H^h = t X) z n [a n a n+1 -l] (lb) 

4 n=-N 

3 

a ■ a — & ' cr 

a=l 

where erf denotes the Pauli matrices operating in quantum spaces V{\ % — — N, . . . , N attached 
to a one- dimensional lattice with 2N sites. 

We choose in ([!]) the ferromagnetic sign of the coupling constant (J < 0) and assume the 
inhomogeneous piece H in h to be a small perturbation of the homogeneous part Hh om , that is, 
we stipulate for the locally varying couplings 

N < \ J\ ■ 

The homogeneous XXX-chain is, as a prototype of an integrable model, one of the most 
thoroughly studied one-dimensional spin models. 

A mathematically rigorous analysis of the model has been provided by Babbitt and Thomas 
|]. For a treatment of the XXX-model in the framework of the algebraic Bethe Ansatz (ABA) 
cf. references [@-||. It is found that the complete spectrum of the model is formed by quasi- 
particles, here called magnons, and bound states of magnons, the so-called string states. The 
integrability of the model implies that the interaction between magnons and strings is of a 
particularly simple structure. It is characterized by the following features 0: 

• multi-particle scattering factorizes into two-particle amplitudes 

• the string states are absolutely stable bound states. 



It follows from these properties that neither genuine multi-particle scattering takes place (with 
a non-trivial reshuffling of the particle momenta) nor does a break-up of the bound states 
occur. 
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We will make use of the technique of the ABA @-||. A first step in this direction is to 
embed the Heisenberg spin model (la) into a family of vertex models. The latter models are 
defined through a monodromy matrix T(A) depending on a spectral parameter A 



T(A) = Ljv(A) • • • L-at(A) 



with the local "Lax-operators" L n given by 

1 



L n (X) = - [2i\l ® l n + a ® <? n ] 



(2) 



(3) 



The unit operator 1 and the Pauli matrices c" act in an auxiliary two-dimensional space, while 
I n and a n act in the quantum space V n . The spin chain model (la) emerges as the logarithmic 
derivative of the vertex model monodromy matrix 



H 



1 d \n(tr T(X)) 

2 dX 



A=0 



2N J n 

n=— N 



2] 



with tr denoting the trace with respect to the auxiliary space. 

The integrability of the vertex models and therewith also the integrability of the XXX-model 
is based on the fact that there is a c-number matrix R = R(X — /£), s.t. the Yang-Baxter- 
Faddeev-Zamolodchikov (YBFZ) relation 



R(X - fi) Li(X) ® Liiji) = Li{n) ® Li(X) R(X - 

is satisfied. R = R(X — /i) is in the case at hand given by 

/ /(^,A) 

R(X -n) = 



(4) 



V 










#(M,A) 1 

1 g(p,X) 










(5) 



fM J 



with 



fM = 1 + 



ic 



fi — X 



and 



ic 



/i — A 



The parameter c is set to unity in the following, meaning that the spectral parameter is taken 
as a dimensionless entity. The local relation (£|) induces the global relation 



R{X - fi) (T(A) ® T{fx)) = (T(n) ® T(A)) i?(A - ^) 



(6) 



which might be considered here as hallmark of integrability. 

The YBFZ-relations can be maintained in certain inhomogeneous generalizations of the above 
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models. One possibility is to choose different representations for different sites of the lattice. 
Such cases have been analyzed in references 0-0. A conceptually simpler possibility, no- 
ticed in [f|, consists in attaching local parameters Zi to the local Lax-operator, that is one 
substitutes Lj(A) by Lj(A — Zj) and obtains then also a modified monodromy matrix 

—N 

T(A, {zi}) = T(X; z„ N ...,z N )=H U(X - z t ) . (7) 

i=N 

One may easily see that the YBFZ relations remain intact, 

R(X - (j) (T(A, { Zl }) ® T(fj,, { Zl })) = (T(jm, { Zi }) ® T(A, {*})) R(X - (j). (8) 

It should be noted that equation (H) only holds in general if one specifies for T(X, {zi}) and 
T(/i, { Zi }) the same distribution of local parameters {z_ N , • • • , z^}. The physics of the model 
on the other hand appears to be invariant under permutations of the parameters. This is a 
consequence of the fact that the BA equations, to be mentioned shortly in the following sec- 
tion, which provide the spectrum of the eigenstates, are insensitive to these permutations. It 
is true (modulo some inessential caveats) that the order of the different representations along 
the lattice, as mentioned above, is for the same reason irrelevant. We therefore believe that 
genuine effects of inhomogenities can only be realized outside the class of integrable models. 
It is an easy undertaking to arrive from the inhomogeneous vertex model (|7|) at the inhomo- 
geneous spin chain (lb). Let us make for this purpose the specifications 

A — ► e A Zj — > ezj . (9) 

The inhomogeneous Heisenberg magnet (lb) is recovered as the logarithmic derivative of the 
vertex model 



H = -J~ln{tr Q T{eX,{eZi})) 
2 de 



J 



N 



e=0 2 



o£( A -*n)l 



n=l 



J X 

= I E( A -^)[^n+i-l]=:AH -H 1 . (10) 

4 71=1 

By taking the derivative with respect to a parameter e which parametrizes different distri- 
butions, one expects to leave the realm of integrability. We will confirm this expectation by 
evaluating non-vanishing irreducible multi-particle scattering amplitudes which are supposed 
to vanish identically in integrable models. 

We will restrict our considerations to the scattering of elementary magnons (the incorporation 
of string states is technically definitely much more cumbersome). A simplifying aspect of the 
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problem can be found in that the inhomogeneous perturbation respects the same global SU (2) 
invariance as the homogeneous term. This implies magnon number conservation. 



The plan of the paper is as follows: in the subsequent section we recall some ingredients 
of the ABA, discuss the thermodynamic limit and introduce the so-called multisite formalism, 
taken from [10]. In section III we evaluate in first order perturbation theory multi-particle 
amplitudes at small momenta. Of crucial importance to achieve this goal will be the repre- 
sentation of form factors as deduced in |L(J. The technical tools to be applied in our analysis 
are approximately the same as those used in [11] for perturbative calculations in antiferromag- 
netic environments (while of course our calculations are simpler and more simple-minded). 
The concluding section is devoted to a qualitative discussion of other kinematical regions of 
multi-particle scattering and a summary. In an appendix we report a perturbative calculation 
of the spectrum of low lying states which is conceptually not related to the theme of the bulk 
of the paper but on a technical level rather similiar to the calculations in section III. 



2 Basics of Bethe Ansatz 



2.1 Algebraic Bethe Ansatz 



We collect here for the sake of self-consistency of the paper some of the basic aspects of the 
ABA \ 

Let the monodromy matrix (|2]) be parameterized as 



T(A) 



A(X) B(X) 
C(A) D(X) 



(11) 



One deduces from the YBFZ relations @ the commutators of the operators A(X), . . . , D(X). 
Of these 16 relations we only list the following 



B(X), B(fi 
B(X), C(fi) 
A^)B(X) 
D^)B(X) 
C(X)D(fi) 







C(A), C( M ) 



g(X,fi)(D(X)A(v)-D(n)A(X)) 
f^,X)B(X)A(^ + g(X,fx)B^)A(X) 
/(A, fi)B(X)D(ii) + g(n, X)B(v)D(X) 
f(\,fi)D(jj,)C(\) + g(ji,\)D(\)C(ji) 



1 For a more thorough introduction see B 
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C(\)A(n) = f(n,\)A((i)C(\) + g(\,ti)A(\)C(fi) . 



(12) 



Let |0) denote the state of highest weight with respect to the tensorproduct of SU{2) repre- 
sentations in the configuration space V = n® Vi- This state is annihilated by the operators 
C(A) and is an eigenstate of the trace of the transfermatrix t(A) 

r(A) = tr T(X) = A(X) + D(X) 
C(A)|0) = 

(A(X) + D(X))\0) = 



lX+ 2) + { lX -2 



10} 



The ABA renders a representation of the eigenstates of the transfermatrix in terms of the 
operators B(X) - being the hermitian conjugates of the operators C(A) - which act as creation 
operators of quasi-particles (magnons) on the highest weight state |0). 
Introducing the notation 

|$(A 1 ,...,A,))= l[B(X n )\0) (13) 



n=l 



one arrives - exploiting the commutation relations 



at 



(A(X) + D(X)) fl B(An)|0> = A(A;Ax, ...,A,) {[ B(X n )\0) + £ k k {\ u ...,A,)S(A) II B(X,)\0) 

n=l n=l k=l j^k 

(14) 

with 



l\ m 'A j-X-i l\ 2N < X 3 - X + i 



(15) 



and 



A fc (Ai, ...,A/) = g(X k , A) 



2N I 



iX k + 



ly | r a, - a/, 



d 5 



A-; — Xi. — i 



2 J J=1 Aj — Afc 



- NA 



j n Ai ~ Afc+i 

2/ J=1 Aj — Afc 



(16) 



The second bunch of terms on the right handside comes from the exchange term in the com- 
mutation relations. The spectral parameters have to be specified s.t. these terms vanish. This 
gives rise to the Bethe Ansatz equations 



' • \ , i \ 2N i \ \ 



iX k - |_ 



n : 

i=l \ - *k - 1 



(17) 
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|$(Ai, ...,Aj)) is an eigenstate of the transfermatrix if the last equations are satisfied. 
The eigenstates can be classified with respect to the eigenvalue of S 3 : 

IN 

s 3 mx 1 ,...,x l )} = (—-i)mx l ,...,x l )} . (is) 

The dual wave functions are given by 

i N 

mm = <oi n ^uo = (-iy<oi n . m 

3=1 J=l 



2.2 Thermodynamical limit 

We introduce suitably normalized operators in order to deal with finite norms of states and 
finite eigenvalues in the thermodynamical (TD) limit ||12||p| : 

i(A) = a-\X)A(X) , D{X) = cT\X)D(X) 
B(A) = a _ 3(A)rf-s(A)B(A) , C(X) = <r*(\)<r*(\)C(\) (20) 

with a(A) = a(X) 2N = (iX + \) 2N the eigenvalue of the operator A(X) and d(X) = 5(X) 2N = 
(iX — \) 2N the eigenvalue of D(X) with respect to |0). Thus the operators A(X) and D(X) have 
correspondingly the eigenvalue 1. Exploiting the relations f] 

lim 9 (A-,)expiM^M = _ rf(A _,) (21) 



and 

1 1 



± in 5(x) 



x x ± ie 

which hold in the sense of generalized functions ||13| - not pointwise - we obtain in the TD 
limit the simplified relations 

A(X)B(fi) = f4X^)B(fi)A(X) 

C(X)D{ft) = f.(X,^)D{fi)C(X) (22) 

with /_(A,/i) = 1 + 3^. 

One notes, comparing with equation (p!2|) , that the exchange terms have dropped out. The 

2 Our prescription differs from the one given in jl^] in an inessential way. 

3 Terms of the form - have to be evaluated according to the principle-value prescription. 
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BAE can therefore be disregarded in the TD limit (as long as one restricts the attention to 
the sector of elementary magnons). Normalized asymptotic scattering states are generated by 
acting with creation operators Z(X) = B(X)A~ 1 (X) on the vacuum (the highest weight state) 
and are annihilated by operators Z'(X) = — D _1 (A)C'(A) pi. The action of the operators A -1 
and D~ l are easily deduced from the relations ( |22|) and from the fact that the vacuum is an 
eigenstate with unit eigenvalue of ^4(A) and D(X) and therefore also of A -1 and D~ x . 
An incoming scattering state is given by 

Z(A x )...Z(A n )|0) (23) 

if the rapidities are ordered in such a way that Ai < . . . < A n , and represents an outgoing state 
for Ai > . . . > A n . 

To relate the incoming to the outgoing states, use has to be made of the relation 

Z{X)Z{n) = S{\, l i)Z{ii)Z{\) (24) 

with S(X, fi) = j^yj the two-body S-matrix. It is easily seen that the n-magnon S-matrix is 
given as a product of 2-magnon S-matrices. 

The wavefunctions (|23| ) are normalized to delta functions with a unit prefactor. 



2.3 Multisite formalism 



To evaluate scattering amplitudes in the Born approximation, we have to determine formfactors 



lTJ of the type 

<oin^(^n^)io> (25) 

i j 

where the operator O is given by O = J2n=-N z n&n&n+i = H n ^n,n+i- So we are lead to 
consider matrixelements of the form 

(oi n z\xf)o n>n+1 n z(xf)\o) = n rw, a?) r w , xfm n c(xf)o n>n+1 n b(\* m 

i j i>j i j 

(26) 

where the latter identity is a straightforward consequence of the definition of Z and the com- 
mutation relation ([22]). 

The basic strategy for the determination of the r.h.s. of (p6|) will consist in decomposing the 
monodromy matrix into parts as follows: 

T(A) = T(3|A)T(2|A)T(1|A) (27) 
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T(1|A) = L n _ 1 (A)---L_ iV (A) 
T(2|A) = L n+1 (A)L i (A) 
r(3|A) = L JV (A)--- J L n+2 (A). 

The sub-monodromy matrices may be parametrized as T(A) above 

r <^> - ( c'M $8 ) 

with(j=l,2,3). The product in (f27j) is meant to be ordinary matrix multiplication of 2 x 2 
matrices. The T(j\X) fullfill the global YBFZ commutation relation seperately, acting on the 
vector space with highest weight |0)j. The highest weight state of the total space is given as 
a tensorproduct 

|0> = |0>3®|0>2®|0>i . (28) 

Using the commutation relations, the operators Aj, Dj, which appear if the B(X) are expressed 
through operators of the subspaces, can be commuted through to the vacuum. This yields the 
so-called multisite formula || 

n*(Af)io>= e ri n ri k b 

3=1 {x BI }u{x BII }u{x BIII }={\ B } j B ei k B eii i B ein 



B 3 (\? B HI )\0) 3 ® B 2 (A£ J )|0>2 ® B 1 (\f;)\0) 1 (29) 



with 



Kb = a 2 (Af s ')<i 1 (Ag')a 3 (A«')<i 1 (Af,''')a 3 (A£')4(AS'") 

x/(Af/,Af'')/(Af/,Af/»)/(Af B ",Af™). 

The summation in (p9|) is with respect to the partition of the set of all Bethe parameters {Xj} 
in 3 disjunct subsets {X BI }, {X BI1 } and {X BHI } with 

card{X BI } = If, card{X BI1 } = l B , card{X BIH } = l B . 

A similar representation can be derived for the dual vector (0| IIj=i C(A?) 

I l l l 2 l 3 

(om^f)= e u u u Kc 

j=l {\ci}\j{\cii}\j{\cni}={\c} j c<E i k c eli laein 

(0\sC 3 (Xg ni ) ® (0| 2 C 2 (Agf ) ® (0|iCi(A£f) (30) 
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with 



Kc = d 3 (A^)o 1 (Ag / )d 8 (Ag)a 1 (A£ m )ifc(Ag J )a a (A£ m ) 



x/(Ag J ,Ag)/(A^,A^)/(Ag" I A£") 
Inserting (EDI) and (BO) into (EH) we obtain 



/ 



(0\UC(\f)O n>n+1 l[B(\Z)\0)= Ell II «j(A£)«MA?)a*(Aj)d.,(A*) 

i=l k=l I,II,III I,II,III I<J<K<III 

x f(Xj, X C K ) /(A?, Af ) 4 /} ({Af }, {Af }) 4 m) ({Af /7 }, {Af /7 }) (0|C 2 (Af 7 )O nin+1 B(Af 7 )iO) 

(31) 

with <Sf.({Af }, {Af }) =j (0| Hi C(Af) Ili S(Af )|0)j being the scalar product in the j-th space. 
The cardinality of the partition sets {Af } and {Af } are equal to each other. Matrixelements 
with card{X^} ^ card{X^} vanish. 

Taking into account the normalization of the operators B and C relative to B and C we arrive 
at the following expression 

<0|ri<7(A?)O nin+1 nS(Af)|0) = 

3=1 k=l 

E II (OlC^A^O^BaCAf^lO) \{ /(Af, Af) /(Af , Af ) 

1,11,1111,11,111 I<J<K<III 



x a5(Af)-^a5(Af // )-^ i a5(Af)-^a(5(Af // )- i ?a(5(Af / )- 1 a(5(A^) 



N 3 



X 



r(Xf 

KAf) 



'(Afj 



iiu 



2 ^ 



r (Af/) 



r(Xfu) 

where aS(X) = a(A)<5(A) and r(A) = 



^ ) ({A?},{Af})4 m) ({Af // },{Af // }) 



(32) 



5(A)- 



One has the following recursion relation for the scalarproducts Jl0|: 



/ 



/ 



S({A?}, {Af }) = a(A? ) E d(A^)^(Af , Af ) [] 0(A? , Af ) II s(A?, Af )S l - 1 (a 1 (X),d 1 (X)) + 



n=l 



i 



k^n 



d(Af ) E «(Af )^(Af , Af ) n ^(Af, Af Af ) n <?(Af , K )^-i(a 2 (A), d 2 (A)) 

n=l jf^l fc^n 
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with 01(A) = a(X)h(X, Af ), a 2 (A) = a(X)h(X, Af ) and di(A) = d(A)/i(Af , A), 

4(A) = d(X)h(X%, A), while fr(A, /u) = j[£4 = 1 + We have quoted here on the r.h.s. the 

functional dependence of the scalar products on the vacuum eigenvalues which have changed 

going from the l.h.s. to the r.h.s. from a(A) to a(A) and d(X) to d(X) respectively, which makes 

the solution of the recursion relation difficult in general. The two-term recursion relation 

simplifies in the TD limit, if we concentrate on the irreducible part of the amplitude. 

We get for the normalized scalarproduct in the limit N — > oo: 



iim nn 




a5(X°)-^a5(Xf)-^Sl?({X?},{Xf}) 



h 



-7T e s(xi 

7ll = l 



^(A^A^n^Af) n gK 1 ,^ 
-A^n^Af) n 9 (*%,w 



ci 



h h 

iim n n 

N^oo J -- L A - L 
00 ci^i Bl^m 



r(Xf 



N 3 



a5(A / c )-^a5(Af)-^4 /) («(A)/ i (A,Af i / ),rf(A)/i(Af / )) 

(33) 



where we have used relation fl2T|). One should note that the second term on the r.h.s. only 
contributes - due to the appearance of a delta function - to scattering processes where at least 
one magnon goes over unscattered from the incoming to the outgoing state. The irreducible 
scattering amplitude however refers by definition to that part of the amplitude from which 
all energy conserving subprocesses have been subtracted. It means that the restriction to the 
irreducible amplitudes effectively implies that only the first term of the recursion relation (|33f) 
has to be taken into account for the evaluation of . The ensuing one-term recursion relation 
is easily solved with the result 



h h 

iim nn 

N ^CIBI 
h 



KAf 
r(A?) 



aSiXfy^aSiXfr^S^iW}, {X?}) 



E E ••• E 



rii=l "2 =1 
™2^ n l 



n ^p,Af) n ^(Af/, xg) . . . n 

j>i=l fci^ni ki^ni,...,™,^ 



ri [MA^,A^)MAS,An]n 



k>i=l 



k=l 



B\ 



r(Xf 
KAF) 



iV3 
2 



<?(A 



CI \BI\ 
k ! A n h ) 



(34) 
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with the subscript "irr" indicating the restriction to that part of the scalarproduct that con- 
tributes finally to the irreducible amplitude. The extension of the r.h.s. of equation (|34j) by 



factors 1 = jfxfy enables us to represent it as a determinant multiplied by some overall factor: 



h h 



^nn 



iV— »oo 



CI BI 



B\ 



r(Xf 



a5(Xf )-^a5(Xf }, {A? }) 



n ^f 1 , Af ) n 9{>? , Af o n ma? 7 , Af ) n 

j>i j>i 



h3 



r(X 



BP 



-, ^3 



r(X< 



det h M {I) (35) 



with M. 



(/) _ gW,\?) 



A similar relation holds for the part denoted with /// (here only the second term in the 
recursion relation of the scalarproduct contributes): 



iim n n 

/V^oo ±A - ±L 

cm Bin 



r(X 



in) 



r(X 



in) 



+i 



u8(X°)-^a5(Xf)-%cals\l n \{X m },{\fu})i 



with M 



n ^(Af 7/ , a? m ) 5 (Af Af") ft fc(Af j ", Af ") ft 



r(A 



KAf") 



+1 



det h M {III) (36) 



Ma 77777 7a7 777 ) • 



Inserting (p5|) and (|36| ) into ( pl| ) one obtains 
(0|nc7(Af O n , n+1 {[B(X*)\0) irr = 



i=i 

n 



it=i 



r(Af) 
r(Af) 



n^(Af AfMAf.Af) £ (-l) [Pc]+[PB] (Oin^(Ajj)a,n + iIl^(Afj)|0) 

7 >i 1,1 1,1 1 1 II II 



h 



II MA|, A?) MAf , Af) J] MAf , Af ) n KX*» Af ") 

KJ<K<IH 



''J 



'•J 



X 



det (R-^X^M^iXf, Xf)R(Xf)) det (R(X? n )M {m) (Xf n , Xf^R'^X 



B 
III 



(37) 



with R(X)ij = r(Xi)Sij. While deriving this result we used /(A,/j) = g(X, fi)h(X, fj) and the 

antisymmetry of the g's. 

[Pb] stands for the parity of the permutation 



P B : {Af,} U {Af } U {Xf n } - {A*} 



11 



while [Pc] stands for the parity of the permutation 

P C : {A&} U {A? } U {xf n } - {X c } (38) 
with the enumeration in each subset according to the original one. 

It is possible in principle to write the result in a more compact way, namely as the determinant 



of the sum of three matrices |L4] . As it is not useful for our purpose we will not pursue this 
line of reasoning. 



3 Low-energy limit 



To start with let us make the simplifications which are due to the special form of the pertur- 
bation. The matrixelement 

2 (0| II C 2 (Xf)a n a n - +1 J] B 2 (Xf )|0} 2 (39) 

i i 

appearing in equation (37) is to be evaluated with respect to the two-site highest weight state 



|0)2. There can therefore at most two operators B and C show up in ( |39"D (applying two 
operators B 2 to 1 0) 2 one reaches the state of lowest weight of the two-site vector space). Since 
we are anyhow restricting our attention to scattering events in non-forward direction we may 
evaluate instead of ([39]) the matrixelement 

2 (0| I[C 2 (X? ) (a n a n - +1 - I n ■ l n+1 ) \{B 2 {Xf )|0> 2 (40) 

i i 

with l n the identity in V n (the addition of l n ■ l n+ i gives only a contribution to the amplitude 
in forward direction). But fl40|) vanishes on the state of highest weight (no operators B 2 
and C 2 ) and on the state of lowest weight (two operators B 2 and C 2 ). We are left with the 
matrixelement with one operator B 2 and C 2 , which is straightforwardly calculated 

2{0\C 2 (Xf I ) (a n a n - +1 - l n ■ I n+1 ) B 2 (Xfj)\0) 2 = 2 1 1 . (41) 

Taking the normalization and the last result into account we obtain for the transition amplitude 
the representation 

E<oiri^(A, c )(v + i nz(A?)io>= 

n 3=1 fc=i 
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r(A; 



r ' ! 2 A g(jf»Af) g(X?,Xf) [Pc]+[Ps] 



/(Af , Af ) /(A?, Af ) , ^ 7 y «5(Af,) «5(Af 2 ; 



n Kx c K , x c j) k(Xj, x b k ) n Kxf 1 , Af ) n ^ m , x?») 

I<J<K<III i,j i,j 

x det (R-\Xf)M^(Xf, Xf)R(Xf)) det (R(X? a )M^(X? u , Af^iT 1 ^)) (42) 

The slash on the sum over the partitions is supposed to indicate that only partitions with 
exactly one representative present in the subset labeled by II are to be taken. 

We are now prepared to examine the behaviour of irreducible scattering amplitudes at low 
momenta (Aj ~ Pi for small momentum) with two or more than two magnons involved (the 
one particle amplitude will be quoted below for the sake of completeness). An obvious method 
to get a handle on formula fl42| ) consists in a systematic expansion in powers of momenta, as 
far as they appear in functions h and keeping at the same time the functions g unexpanded. 
The leading term is obtained by putting h consistently to one at all places where it appears in 
(|2|). This yields fot A Q e {A c } 

£<O|rizt(A?)0 nin+1 I3Z(Af)|0>« 

n j=l k=l 

£32^n ^'^ff ^ Y.'^r^detgiXfA^detgiX^A?!!) • 

j>i J( A j ' A i I A A i ' A 3 ) I, II, III 

(43) 

The prefactors 0j->* jt^jP) ma y & l so be put equal to one in leading order by noting that 
yjy^y ~ 1 + O (X — fj,). The sum over the partitions I and III in (|43| ) renders a vanishing result 
as one infers from the Laplace formula for the determinant of a sum of matrices ||: 

det{A + B)= (- l ) [PL]+[Pc]detA PL,Pc detB PL,Pc ( 44 ) 

Pl,Pc 

where Pl is the partition of rows in subsets of rows of A and B, while Pq is analogous the 
partition of columns, and the fact that g is odd 

g(Xf,Xf)+g(Xf,Xf) = 0. 

For the next order of the expansion in powers of momenta we obtain in a straightforward 
manner the following result: 



(0| II Zt(A?) af + i n Z(X*)\0) « 2t {[ fr\\f, Xf^Xf, Xf) J] 

j=l k=l 3>i i 



KAf) 



KAf) 

16 4 £ A Q £'(-l) [ft]+[PoI (fe - k + Qe 1 / 11 ) detg(Xf, Xf) detg(Xf Il7 X? n ) (45) 



x 

c 1,11,111 
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with 



j,iii 



+1; K e / 
-1; A a G /// 
0; A a G 77 



If A a G {A B } we get the same result up to an overall minus sign. 

The result can be simplified further, using the following chain of identities: 



£ (_i)[^M^](Z s - / x ) detg(Xf, Xf) detg(Xf H , X 



in) 



ijii 



d 



7T E (" 1) [Psl+[Pcl 9(xXf, xXf) det g{x- l XU z -1 A?„) U 

^ 7,777 



Q (\ \ h+h 



= — det (g(xX c , xX B ) + g{x~ 1 \ B , x~ x X 

detg{X c ,X B ) 



c 



x=l 



dx \x 
; for h + l 3 > 1 



where we used again the Laplace formula and the antisymmetry of the <?'s. 
The remaining term of the first-order Taylor expansion is 

<0| fl Z^(Xf) Zi a- af +l n Z(Af )|0> = {[ h~\Xf , Xf)h~\X B , Xf) J] 



r(A 



x 16 -6iziJ2 A « E'(- 1 ) [Psl+[Pcl (4' m )^^(Af,Af)de^(Af J/ ,A^). 

C 7,77,777 

The sum over the first and the third partition can be combined to a determinant of the sum 
of two matrices 



E (-l) [Psl+[Pcl (4' m ) 0(Af , Af ) det »(Af„, A 

7„777 

= det (g(X c ,X B )+g a (X B ,X 



in) 



The matrix g a (X B , A c ) differs from ^(A- 6 , A c ) in that the ct-th row is multiplied by (—1). The 
sum of the two matrices is thus a matrix with only one row of non-vanishing entries. The 
determinant is again zero except for h + h = 1- We keep as net result that the first-order term 
of the Taylor expansion gives a non-vanishing contribution only for the one- and two-particle 
amplitude. 

We turn now to the second-order contribution. The computation is tedious, but it proceeds 
otherwise along the same lines as the first-order calculation. We thus only quote the result 



14 



(omitting the prefactors): 

(0|nzt(Af)^af + i II^(Af)|0> 



3=1 k=l 
1 

2 



16 \ E E'(-l) [Psl+[Pcl [*£' C ('i " *3 - ^)ief c (h 6s 1 / 11 ) 



a/3 1,11,111 

~i.ni 



+M34 e~ J ' 7 " + h + h) detg(Xf, Xf) detg(Xf n , Xf n ) (46) 



with 



and 



b,c _^ J +i; x a <e b 

-\ — 1; A Q G C 



_r, m= f +1; A in/,/// 
1 0; otherwise 



e 



When Zi H- Z 3 > 2 this term vanishes as can be shown by generalizing the considerations 
used with the first-order calculation: 

• terms proportional to (/ x + l 3 + const) vanish by the same argument as used in the 
first-order calculation (even for l 1 + l 3 — 2) 

• terms proportional to — Z3) 2 vanish when regarded as the second derivative with respect 
to x at x — 1 (this term renders for l± + Z 3 = 2 the only nonvanishing contribution) 

• terms proportional to e^ ,m (/i — l s ) give a matrix with a prefactor — 1 3 (| + 
of which the derivative with respect to x at x — 1 vanishes 

• terms proportional to e^' m e^' m give a matrix with at most 2 columns or rows or one 
column and one row not zero after applying the Laplace formula 

We finish this section by quoting the leading terms of the transition amplitudes at small 
momenta with the explicit expressions for I < 3 in lowest order: 



I = 1 



Y,z n (0\Z\X)a n a n -' +1 Z( f x)\0) = -2^>n 



N±-N 3 

r(/i)_ 



A 2 + i M 2 + 1 



(47) 
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1 — 2, cf. equation ([45] 



J2z n (0\Z^X^Z^X^)a n an +1 Z(X^Z(X^\0) = 

n 

-64 E^ (EA C -E^) 2 det 2 (^1 



x B 



(4* 



I = 3, cf. equation ([46] 



" i=l i=l 

-128 E ^ (e A ° - E AB ) 3 rfrf 3 (r^ 

n \ c B / ^ 



X B 



(49) 



with deti ( juz^b ) denoting the Cauchy determinant of a I x I matrix 
det, ( -^—) = (-1)^ 



1 \ , 1 ni<,-(Ai - Xj) Ui<j(fH - Hj) 



A couple of remarks may be in order: 

1. The one-particle amplitude quoted above is in fact the full Born term (not the leading 
piece at small momentum). 



2. We have left out in ((48]) and (g|) factors Hi [7M} 2 ■ The latter g ive rise to a 
Fourier transform of the distribution of coupling constants. We have effectively evaluated the 
distribution at zero momentum. 

3. To apply the above expressions to physical processes of magnon scattering one has to 
restrict the respective expressions to the energy shell, given by Ya A f 2 — Yi ■ 



4 Conclusion 



The main result of this paper are the formulas ( 4"8|) and fl4*9|) for two- and three-magnon 
scattering at small momenta. An obvious generalization to n-particle scattering may be con- 
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jectured: 

(0|nzt(Af)a>f + i f[Z(X?)\Q)K-162 l ^X c -^X B ^ ^(jc^s) • ( 5 °) 
We are not able to prove this conjecture so far. 



There are other kinematical regions besides the one of low momenta for which simple and 
reliable estimates can be made. If all momenta and all differences of momenta become large, 
the n-particle transition amplitude decreases with denoting; a common scale of all 

momenta - as can be inferred from an inspection of equation ([37]). An interesting kinematical 
region - also accessible to a rather detailed analytical description - is given by the setting 

|Af-Af|«l, |Af-Af|<l, |Af - Af| > 1, ,\X B \ ~ \X C \ ~p> 1. 

This situation is realized if a bunch of particles travelling approximatively with the same 
velocity is collectively scattered backwards at the inhomogeneity. 
The piece of fl37| ) supplying the p dependence in this case is given by 

a (0|C a (Ag) £ 2 (Af 7 )|0> 2 ft KXf J , Af ) ft h(X Bm , Xf n ) 

i,j id 



det (M^\Xf, Af )) det (M^{X? ni \f n ) 



for which one easily calculates the scaling behaviour p~( n + 3 ). 

To arrive at this conclusion it is essential to view the determinants in the above formula as 
derivatives of Cauchy determinants: 

detMij * detl JxTx-f ~ W ' ' ' M detl Jx^x~) 



A completely open problem within our approach is the treatment of string states. The 
determination of break-up amplitudes for string states seems to us a particularly challenging 
problem. 
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A Appendix 



In this appendix we use the result fl31|) to determine the shift of energy eigenvalues caused by 
H inh = \ Eli z n [a n a n+1 - 1] p|, 0. 

The lowest excitation is generated by flipping one spin (7 = 1). The solution of the Bethe 
Ansatz equation is in this case 



\ 1 f P° 
A = - cot — ; 

2 2 ' 



27rk 



k= 1,...,N . 



(51) 



Taking into account parity degeneracy the first-order correction to the energy E^°'(X) = ^ x2 \.i 
is found to be 

E 0) 



V(A,A) ±2 _ |V(A,-A)| 



(0|C(A)S(A)|0) " (0|C(A)B(A)|0) 



with V(/zi,// 2 ) = — \ Ej=i - 2 i(0|C'( / ui)((Tj(Tj + i — I)B(yU 2 )|0) which leads to 



E^\\) = E^{\) 



N 



N 



£ "-0 ex P (-2zp (A)(i-fc)) 

\ j,k=l 



(52) 



This shows that the energy correction depends in first-order both on the mean-value of the cou- 
plings z = jj Z)j=i z j an d on the Fouriertransform of the distribution Y.f=i z j ex P(=t2ip j) 

The second-order corrections can be obtained from the secular equation 

^nrri^mn 



4 2) (A) = E 



- P (o) P (o) 

m^n -Cm C/m 

where the matrixelements are taken with respect to the corrected wave-function in the zeroth 
approximation 

V nm = (cS 0) ^(A)±4 0) ^(-A)| J f/ 1 |$(^)) 

= ci°h;(A^)±4 0) *v(-A,/i) 

with V(A,/x) defined as in (^) and c*- -* being the following expressions 



c (0) 



c {0) 

c 2 



V(A, -A) 
\]2|V(A,-A)| 

V(-A,A) 
\\2\V(-\,\)\ 
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Inserting the explicit formulas yields 

1 



4 ±t, \ 2 -^ 2 N 2Z -t 

/i^±A 3,k 

j-k 



^l Z j Z k 



1/jL 



j-k 



± 



K± is a quotient of Fouriertransforms 



iX + \ 
Jx^\ 



ifJ>+%J 



'iX 



j-k 



j-k' 



iX 



' iX~l 

A + i-, 



iX + 



j+k- 



^=1 Z J I a T i 



iV 



\l ^i=i ^ ^Slf 



2.7 



2.7 



The sum over /x can be transformed for N — > oo into a principle-value integral 



A 2 - /i 2 



Ci-fc-i) 



j\ (fc-i-i) 



""2 



2cos[p(A)(i-A)] 



" + 2 



ir+e -ip(A)(i+fe) + K _ e ip(A)(i+fe)l I _ 



This principle-value integral can be evaluated by deforming the integration contour into the 
complex plane, closing it at infinity, which is possible as the integrand vanishes as r -4 at 
infinity. Thus only the pole structure of the integral matters. 

There are poles at A = ±/i for all values of j, k, at /j = | for j > k in the first term and for 
j < k in the second term and at /i = — | with j, k dependence of the first and second term 
interchanged. It is convenient to split the sum over j, k into three parts: 

E=E+E+E • 

j,k j=k j>k j<k 

For each sum the contour can be deformed in such a way that the integrand only contains 
poles of first order, for which the residues are easily calculated. 



The result of the integration is 

4 2) (A) 



-£°(A)<^E|[2±Mi = *)] 



7 E ^ sin b(A) (j -k)}[2 cos b(A) (j -k)]±f(X,j> k)} \ (53) 

j>k 



with f(X,j,k) = \K+e- ip ( x W + Q + K„e ip ^ + ^ 
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For the second lowest excitation (two magnons) the computation is more involved, but 
still elementary, so we only give the result for the first-order correction to the energy of the 
two-magnon state: 



£ (1) (/^ 2 ) 



N 



£ (0) (/^A^E^±2 



■ JV 

E 

.j,k=i 



Zj Zfc 



exp [ - 2ip (fi 1 )(j - k)] + exp [ - 2ip (^ 2 )(j - k)] 



N N 

±2£( )(^)£ (0) ME ^exp[-2ip ( A n)(7-A:)] E ^F? ex P I ~ 2<PbM(j " *)] 



i,fc=i 



AT2 



Furthermore there exist complex solutions of the Bethe Ansatz equations. They describe bound 



states 17 with momentum 



D ip(x) 



x + i 



X — I 



and energy 



E String( X ) ~ x 2 + I 



where x denotes the center of the bound state. 

The first-order correction for the two magnon bound state is 



Estringi x ) ~ ^String( X ) 



with Eqm the energy of the center of the bound state. 



1 



N 



N 



4 E -h? exp(-2ip(x)(j - k)) 
\i,k=i N2 



(54) 
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